New measurement on photon yields from air 

and the application to the energy estimation 

of primary cosmic rays 



'^ '. M. Nagano^, K. Kobayakawa^, N. Sakaki*^ and K. Ando' 

O 

^\i ■ ^Department of Space Communication Engineering, Fukui University of 

w\. Technology, Fukui, 910-8505 Japan 

5 ■ ^Department of Architecture and Civil Engineering, Fukui University of 

Technology, Fukui, 910-8505 Japan 

^RIKEN (The Institute of Physical and Chemical Research), Wako, 351-0198 

Japan 

Kj^ . '^Department of Applied Physics and Chemistry, Fukui University of Technology, 

'sj- ' Fukui, 910-8505 Japan 

o. 

^ ■ Abstract 

o 

1-C ' The air fluorescence technique is used to detect ultra-high energy cosmic rays 

II . (UHECR), and to estimate their energy. Of fundamental importance is the photon 

O I yield due to excitation by electrons, in air of various densities and temperatures. 

'^ ' After our previous report, the experiment has been continued using a ^''Sr /3 source 

d ■ to study the pressure dependence of photon yields for radiation in nitrogen and dry 

air. The photon yields in 15 wave bands between 300 nm and 430 nm have been 

determined. The total photon yield between 300 nm and 406 nm (used in most 



< 

o 



x 



c^ 



^ • experiments) in air excited by a 0.85 MeV electron is 3.81ib0.13 (±13 % system- 

atics) photons per meter at 1013 hPa and 20 °C. The air density and temperature 
dependencies of 15 wave bands are given for application to UHECR observations. 
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1 Introduction 



In order to detect ultrahigh-energy cosniic rays (UHECR), atmospheric fluo- 
rescence light from the trajectory of the extensive air shower may be measured 
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by mirror-photosensor systems (the fluorescence technique) . In this type of ex- 
periment, the photon yield from electrons exciting air of various densities and 
temperatures is the most fundamental information for estimating the primary 
energy of UHECR. We have reported in our previous paper [1] that the pho- 
ton yield between 300 nm and 406 nm in air excited by an electron with a 
mean kinetic energy of 0.85 MeV (^"^Sr j3) is 3.73±0.15(±14% systematics) 
photons per meter at 1000 hPa and 20°C. We measured it with six narrow 
band interference filters, whose central wave lengths were 314.7, 337.7, 356.3, 
380.9, 391.9 and 400.9 nm. The bandwidth of each filter at 50% of the peak 
transmission was about 10 nm. We estimated the photons in unmeasured wave 
band to be 8.8% with help from the values reported by Bunner [2]. We have 
also shown that the photon yields in the 337 and 358nm bands are propor- 
tional to dE/dx, when we compare our measurements with the photon yields 
measured at 300~1000MeV by Kakimoto et al. [3,4]. 

We have made the new measurements with additional 9 filters and have im- 
proved the photon yields reported before. By comparing the measurements 
with filters of overlapping wavelength, the contamination of bands in the tail 
of each filter is estimated and corrected. In this report we provide the photon 
yields between 300nm and 430nm from our own measurements. The air den- 
sity and temperature dependences of each wave band, and the average values 
in the radiative transition from the level v' to the level v" in the first nega- 
tive 1N(0, v"), and in the second positive 2P(0, v"), 2P(1, v") and 2P(2, v") 
systems, are given for application to UHECR observations. 

In our previous report [1], we have shown the pressure dependence of the re- 
laxation rate of the excited level of N2 and N^ in nitrogen gas and in air. 
The radiative lifetimes Tq have been obtained after determining the reference 
pressure p' from the pressure dependence of the photon yields (see Eqs.(5) 
and (10) in [1]). The definitions of tq and p' will be given in Section 3. Since 
there have been many experiments (e.g. [5], [6]) before the 1970s that mea- 
sure To with accuracy much better than ours, we don't include our present 
measurements of Tq in this article. 



2 Experiment 



We chose a photon counting and thin target technique to measure the pressure 
dependence of photon yields (the number of photons produced by electrons 
per meter of travel) from nitrogen and air excited by electrons, following the 
method employed by Kakimoto et al. [4]. Experimental details are described 
in our previous report [1]. 

The central values of the new filters used in the present measurement are 325.0, 
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Fig. 1. Typical example of quantum efficiency of PMT used (a thin solid curve) 
and transmission coefficients of interference filters used in the present experiment 
(dashed curves). Thick solid curves are those used in the previous report [1]. The 
relative intensities of fluorescence bands at 1000 hPa from the nitrogen molecule 
(thin vertical lines) and ion (broad vertical lines) in air from Bunner [2] are shown. 



330.6, 350.2, 372.5, 410.0, 414.0, 418.5 and 430.0nm, and their bandwidths 
at half maximum are about lOnm. A measurement at 392. Onm with a band 
width of 4.35nm is also made. Their transmissions are indicated by dashed 
curves in Fig.l. The solid curves are transmissions of the filters used in the 
previous measurement [1]. By subtracting the contribution of bands in the 
tail of the transmission curves, the photon yields from 15 bands between 300 
and 430nm are determined. Contributions from the 311.7, 313.6, 330.9, 333.9, 
346.9, 350.0, 358.2, 367.2, 371.1, 388.4 and 389.4 nm bands can't be separated 
and are included in other bands. 



The systematic errors of the present experiment are discussed in [1]. Although 
the uncertainty from the contamination of lines in the tails of filter transmis- 
sions is reduced from 5% to 2%, the total systematic error is only reduced from 
13.8% to 13%. This is because the main systematic errors are from uncertain- 
ties in the collection efficiency and quantum efficiency of the photomultiplier 
tubes used in the present experiment. 



3 Photon yield 



As described in Bunner [2] and in our previous report [1], the photon yield 
e per unit length from gas excited by an electron is written as a function of 
pressure p at a constant temperature T in Kelvin: 

p fdE\ ( ^° \ Cp . . 



RT{hv) \dx) \l + ^ 1 + ^ ' 



where R is the specific gas constant (N2 : 296.9 m^s^^K^^ and Air : 287.1 

m^s^^K"^) , — — is the energy loss in eV kg^^ m^, and hv is the photon energy 

da; 
(eV). $° corresponds to the fluorescence efficiency in the absence of collisional 

quenching, p' is the reference pressure where the lifetime of collisional de- 
excitation is equal to the combined lifetime, tq, of the excited state for decay 
to any lower state, and of internal quenching. 

The fluorescence efficiency for the ith band at pressure p, $i(p), is expressed 
by 
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where <l>° is <l>° of the ith band. 

The p' for nitrogen is written in terms of Tq and the cross-section for nitrogen- 
nitrogen collisional de-excitation ann as 

1 4to 

p y/nMnkT 



where M„ is the N2 molecular mass and k is Boltzmann's constant, p' for 
air is related to a"„„ and the cross-section for nitrogen-oxygen collisional de- 
excitation a no as 
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where Mo is the mass of oxygen molecule, and /„ and fo are the fractions of 
nitrogen (0.79) and oxygen (0.21) in air, respectively. 



It should be noted that since p' is proportional to VT, then jo' at T, p^, may 
be expressed by 



^^=^1^' ^^^ 



where ^20 i^ v' ^^ 20°C. 



4 Analysis 

^.i Derivation oj photon yields 



The photon yield per unit length per electron e is determined as the number 
of signal counts A^ divided by the product of the following: the total number 
of electrons /, the length of the fluorescence portion a, the solid angle of the 
PMT Q, the quartz window transmission rj, the filter transmission /, and the 
quantum efficiency QE and the collection efficiency CE of the PMT. 

(6) 



Ixaxi^xr]xfx QE x CE 



In the following analysis we use a value of / appropriate for the main nitrogen 
emission band in each filter band pass. The number of photons from another 
band in a given filter is estimated from measurements in two adjacent filters 
and is subtracted from the observed photons. / is about (4 ~ 5) x 10® from 
about 80 hours in each run. Finally, e from a run in a vacuum is subtracted 
from each e determined above and the corrected e is determined. 



4-2 Two components analysis method 



In order to separate the photons from the IN band system (427.8 nm) and the 
2P band system (427.0 nm), a "two components" analysis has been performed 
in one filter band. The method is described in detail in our previous report 
[1]. Briefiy, we have fitted the observed pressure dependence of eobs(p) with a 
superposition of two bands in one filter by the least square(LS) method. In 
this case the observed photon yield eobs(p) is the sum of the photon yields of 
the main band ei{p) and the sub-band 62 (p), and is written as follows : 



eobs(p) = ei(p) + e2(p) = — ^ + — ^ , (7) 

where x = -, ai = ^ and 02 = 4-- Ci and p\ are the parameters of the main 
band in the filter, and C2 and p'2 are those of the other band. 

We determine a set of four parameters ai, Ci, 02 and C2 in Eq.(7) by the LS 
method. 



5 Results 



5. 1 Nitrogen 



The pressure dependencies of photon yields in nitrogen gas with fourteen filters 
are plotted by solid circles in Fig. 2. In each figure, the main emission band 
in each filter band is listed. The solid curves in the figures are the best fits 
of Eq.(l) with p' listed in the figure. Open circles are the yields of the main 
band in the filter after the subtraction of other bands, which was determined 
by taking into account the transmission coefficients of other bands. The errors 
are relatively large due to the uncertainties of these subtractions. The values 
of p' and C for each band are determined by a fit to Eq. (1) by the LS method 
as described in [1]. 

Photon yields (e) per meter per electron of average energy 0.85MeV in nitrogen 
gas at 1013 hPa and 20°C are determined with p' and C, which are listed in 
the third and the fourth column in Table 1. $° values are also listed in the 
table. The total e between 300 and 406 nm is 21.69±0.55, and that between 
300 and 435 nm is 22.20±0.56. 



5.2 Air 



The pressure dependencies of photon yields in dry air, which is a mixture of 
78.8% nitrogen gas and 21.1% oxygen gas, in fourteen wave bands are indi- 
cated by solid circles in Fig. 3. Open circles are the yields of the main band in 
the filter after the subtraction of other bands, which are estimated from mea- 
surements in neighbouring filter bands, taking into account the transmission 
of those lines. 

The values of e in air at 1013 hPa and 20°C are determined from Eq.(l) with 
the best fitted values of p' and C. These together with $° values are shown in 
Table 2. In estimating e, C and <l>°, the contamination of lines in the tail of 
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Fig. 2. The pressure dependence of e in N2 at 20° C excited by electrons with an 
average energy of 0.85 MeV in 14 filters of different central wave length or different 
wave bands. In each figure, the main contributing line in each wave band is listed. 
Solid circles are the measured values. Open circles are the values after the subtrac- 
tion of contamination of other emission bands. Black solid curves are the best fits 
of Eq.(l) to the experimental points without subtraction. Red ones are the best fits 
to the corrected points (open circles). In the figure labeled 427 and 428 nm, results 
of a two-components analysis are shown. Open circles with a solid line are from the 
IN band system and the dashed line represents the 2P band system. 



the transmission curve are subtracted from each other as mentioned before. 
Therefore some values listed in this table are somewhat reduced from the 
values in Table 9 of [1], where the measured values in each filter at 1000 hPa 



Table 1 

Summary of the measurements of nitrogen gas. See text for detail. 



main 
A(nm) 


e 


P' 


c 


(|)° 


m-i 


hPa 


xlO-V(hPa-m) 


xlO-^ 


316 


2.03 ±0.21 


88.1 


± 7.5 


2.51 ±0.14 


5.07 ±0.28 


329 


0.622±0.063 


121. 


±10. 


0.575±0.033 


1.12 ±0.06 


337 


8.28 ±0.25 


155. 


± 4. 


6.16 ±0.10 


11.7 ±0.2 


354 


0.417±0.044 


70.3 


± 6.4 


0.634±0.035 


1.15 ±0.06 


358 


5.64 ±0.31 


125. 


± 6. 


5.07 ±0.16 


9.07 ±0.29 


376 


0.873±0.059 


82.5 


± 4.7 


1.14 ±0.04 


1.95 ±0.07 


381 


2.09 ±0.25 


128. 


±14. 


1.84 ±0.09 


3.08 ±0.16 


391 


0.419±0.049 


5.46± 0.50 


7.72 ±0.54 


12.6 ±0.9 


394 


0.185±0.078 


39.4 


±12.5 


0.49 ±0.13 


0.79 ±0.22 


400 


0.399±0.036 


62.9 


± 4.8 


0.674±0.033 


1.08 ±0.05 


406 


0.73 ±0.15 


140. 


±25. 


0.597±0.064 


0.94 ±0.10 


414 


0.108±0.029 


111. 


±24. 


0.108±0.017 


0.167±0.027 


420 


0.073±0.028 


34. 


±10. 


0.222±0.050 


0.338±0.076 


427 


0.188±0.113 


232. 


+144. 
- 71. 


0.099±0.038 


0.148±0.057 


428 


0.151±0.031 


5.6 


± 1.1 


2.72 ±0.24 


4.07 ±0.36 


Sum 


21.69 ±0.55 


(300nm~406nm) 


Sum 


22.20 ±0.56 


(300nm~430nm) 



were listed. Total photons per meter per electron between 300 and 406nm is 
3.81±0.13 and between 300 and 428nm is 4.05±0.14. 



6 Discussion on experimental results 



6. 1 Comparison with our previous report 



The present results are obtained after subtracting the contamination of other 
bands in the tail of each filter, by comparing the measurements with another 
filter of overlapping wavelength. In this subtraction we have used the fitted 
curve to the experimental points. The photon yields at 1013 hPa in Table 



Table 2 

Summary of the measurements of air. 



main 
A(nm) 


e 


p' 


C 


(I)° 


m-i 


hPa 


xlO-V(hPa-m) 


xlO-4 


316 


0.549±0.057 


23.0 ± 


1.9 


2.44 ±0.15 


4.80 ±0.29 


329 


0.180±0.026 


40.2 ± 


4.6 


0.465±0.042 


0.880±0.080 


337 


1.021±0.060 


19.2 ± 


0.7 


5.43 ±0.15 


10.01 ±0.27 


354 


0.130±0.022 


30.6 ± 


3.9 


0.437±0.046 


0.769±0.080 


358 


0.799±0.080 


18.1 ± 


1.4 


4.50 ±0.28 


7.82 ±0.48 


376 


0.238±0.036 


34.1 ± 


4.1 


0.722±0.068 


1.20 ±0.11 


381 


0.287±0.050 


19.4 ± 


2.6 


1.51 ±0.17 


2.46 ±0.27 


391 


0.302±0.020 


5.02± 


0.26 


6.04 ±0.25 


9.60 ±0.39 


394 


0.063±0.033 


24.2 ± 


9.4 


0.267±0.093 


0.42 ±0.15 


400 


0.129±0.019 


24.2 ± 


2.8 


0.544±0.053 


0.847±0.082 


406 


0.118±0.019 


12.3 ± 


1.6 


0.972±0.010 


1.49 ±0.15 


414 


0.041±0.009 


19.3 ± 


3.4 


0.217±0.031 


0.327±0.047 


420 


0.042±0.015 


7.3 ± 


1.9 


0.58 ±0.13 


0.86 ±0.20 


427 


0.032±0.023 


70 +60. 
'^- -23. 


0.047±0.021 


0.069±0.031 


428 


0.121±0.022 


3.86± 


0.59 


3.14 ±0.28 


4.57 ±0.41 


Sum 


3.81 ±0.13 


(300nm~406nm) 


Sum 


4.05 ±0.14 


(300nm~430nm) 



1 and 2 are calculated with the best fitted values of p' and C . In contrast, 
in our previous report, the measured values at 1000 hPa are listed and the 
unmeasured values are estimated from the table by Bunner [2]. 

In Fig. 4, the present values of e per meter per electron at 1013 hPa and 20°C 
are compared with the previous results. The closed circles are the present 
results and the open circles are the sum of two bands where they were not 
separated in the previous measurement. The lines are drawn to compare them 
easily. There are slight differences between the new and old results at 329 
nm and 391 nm. In the previous report we estimated the value at 329 nm 
from that listed in Bunner [2], which is relatively low compared with the 
present measurement. Considering the 391 nm value, in the previous report 
we separated the 394nm band from the 391nm band by using a two components 
analysis, while in the present experiment the filter of the narrower bandwidth 
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Fig. 3. The pressure dependence of e in air at 20°C excited by electrons with an 
average energy of 0.85 MeV in fourteen wave bands. In each figure, the main con- 
tributing band in each filter is listed. Solid circles are the measured values. Red 
open circles are the values after the subtraction of contamination of other bands 
described in the text. Black solid curves are the best fits of Eq.(l) to the experi- 
mental points. Red ones are the best fits to open circles. In the plot labeled 427 and 
428 nm, the results of a two-components analysis are shown. Open circles with the 
solid line are from the IN band system and the dashed line represents the 2P band 
system. 

is used to separate them. Though the combined yield at 391nm and 394nm is 
similar, the individual values are somewhat different as seen around 391nm in 
Fig. 4. 
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Fig. 4. Comparison of the present photon yields per meter per electron (closed 
circles) with our previous report (open squares) [1] . Open circles are the sum of two 
bands which were not separated in our previous measurement. 

6.2 Density and temperature dependence of the photon yields 



The photon yields per meter by an electron of energy E can be rewritten as a 
function of the gas density p (in kg m"'^) and the temperature (in Kelvin) to 
apply to the atmosphere [4]. 



(M) 



Ap 



(^) 1 + pBVT ' 

V dx y o.85MeV 



where 



A 



(M) 
V dx y 



85McV 



$° 



hu 



(9) 



and 



Rn 



B = -^\/293 . 
F20 



(10) 
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Table 3 

A and B values in 15 bands 



main 
A(nm) 


Nit 


rogen 


Air 


A 


B 


A 


B 


2t -1 

m^kg 


m^kg-iR-^ 


m^kg-i 


m^kg-iR-^ 


316 


21.8 ±1.2 


0.577±0.049 


20.5 ±1.3 


2.14±0.18 


329 


5.00±0.29 


0.419±0.035 


3.91±0.35 


1.22±0.14 


337 


53.6 ±0.9 


0.328±0.008 


45.6 ±1.2 


2.56±0.10 


354 


5.52±0.31 


0.723±0.066 


3.68±0.39 


1.60±0.21 


358 


44.1 ±1.4 


0.407±0.019 


37.8 ±2.3 


2.72±0.22 


376 


9.95±0.35 


0.616±0.035 


6.07±0.57 


1.44±0.17 


381 


16.0 ±0.8 


0.397±0.042 


12.7 ±1.4 


2.53±0.35 


391 


67.2 ±4.7 


9.31 ±0.86 


50.8 ±2.1 


9.80±0.51 


394 


4.3 ±1.2 


1.29 ±0.41 


2.25±0.78 


2.03±0.79 


400 


5.87±0.28 


0.808±0.061 


4.58±0.44 


2.03±0.23 


406 


5.20±0.56 


0.363±0.064 


8.18±0.82 


3.99±0.52 


414 


0.94±0.15 


0.46 ±0.10 


1.83±0.26 


2.55±0.45 


420 


1.93±0.43 


1.49 ±0.46 


4.9 ±1.1 


6.8 ±1.7 


427 


0.86±0.33 


0.22 ±0.10 


0.40±0.18 


0.68±0.38 


428 


23.7 ±2.1 


9.1 ±1.8 


26.5 ±2.4 


12.7 ±1.9 



where Rair is the specific gas constant for air (see equation 1) . The values of 
A and B are calculated and are listed in Table 3. 

Using the values of A and B in Table 3, the photon yield of a 0.85 MeV electron 
is shown in Fig. 5 as a function of altitude. The US Standard Atmosphere 1976 
[7] is assumed in the left-hand figure. Thin dashed lines are from the 2P band 
system and a thin solid line is from the IN band system. A thick line shows 
the sum from all bands between 300nm and 406nm. The altitude dependence 
from Kakimoto et al. between 300 nm and 400 nm is also shown by a thick 
dotted line. Since the proportion of yield from the 406nm band is 3.2%, there 
is a difference of about 14% (300~400 nm) between ours and Kakimoto et al. 
independent of altitude. This is because the photon yields of Kakimoto et al. 
were underestimated in their unmeasured bands. 

The altitude dependence of photon yields in January and July are compared in 
the right-hand figure for latitudes 30°N and 60°N, assuming the US Standard 
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Fig. 5. The altitude dependence of photon yields per meter for a 0.85 MeV electron 
( between 300 and 406 nm). In the right- and left-hand figures, the US Standard 
Atmosphere 1976 and 1966 are assumed, respectively. In the left-hand figure, thin 
dashed lines are from the 2P band systems and a thin line is from the IN band 
system. A thick solid line is the sum of all bands between 300 and 406 nm. The 
relation from Kakimoto et al. is also shown with a dotted line for wavelengths 
between 300 and 400 nm. The photon yields in January and July are compared in 
the right-hand figure at latitudes of 30°N and 60°N. 

Atmosphere 1966 [8]. The annual value is also shown for 15°N. Monthly mean 
values of temperature and densities from CIRA 1986 (COSPAR International 
Reference Atmosphere) lie approximately between the lines shown for January 
and July at 30°N [9]. It would be necessary to use the several types of alti- 
tude dependence at different latitudes in experiments from space like EUSO 
(Extreme Universe Space Observatory) [10]. 



6.3 Average p' values in different radiative systems 



For the radiative transition from the level v' to the level v" in the IN and 2P 
systems, p' values are expected to be the same for a fixed v' irrespective of 
v" , if $° is much smaller than unity. So in Table 4, we summarize the average 
values of p', B and the ratio of photon yields in nitrogen gas and air at 1013 
hPa. Here weighted averages of four bands (337, 358, 381 and 406nm) for 2P(0, 
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v"), five bands (316, 354, 376, 400 and 427 nm) for 2P(1, v"), two bands (394 
and 420 nm) for 2P(2, v") and two bands (391 and 428 nm) for 1N(0, v") are 
taken for each parameter. 

The present p' values are larger than those listed in Bunner [2], except 2P(0,f") 
in air. The reason is not clear. In the next sub-section we compare the fiuores- 
cence efficiency at pressures of 100 ~ 1000 hPa directly with other experiments. 

Table 4 

Average p' and B values in different radiative systems in nitrogen gas and in air. 



transition 
state 


gas 


P' 


B 


e{N2)/e{air) 
at 1013 hPa 


P' [2] 




hPa at 20° C 


m^kg-iR-^ 


hPa at 27° C 


2P(0,z;") 


N2 


144.7 ±3.1 


0.343±0.013 


7.70±0.42 


120 


air 


18.1 ±0.6 


2.72 ±0.09 


20 


2P(l,z;") 


N2 


74.5 ±2.8 


0.628±0.03 


3.43±0.29 


32.6 


air 


25.6 ±1.4 


1.92 ±0.10 


8.7 


2P(2,?;") 


N2 


36.2 ±8.0 


1.38 ±0.26 


1.94±0.81 


14.5 


air 


7.9 ±1.8 


6.2 ±1.4 


6.1 


m{0,v") 


N2 


5.48±0.46 


9.27 ±0.25 


1.36±0.16 


1.6 


air 


4.83±0.24 


10.2 ±0.5 


1.6 



6.4 Comparison with other experiments 



The present values of ^i{p) given by Eq.(2) in N2 and air at p =800 hPa (~600 
torr) are listed in Table 5, together with those of Davidson and O'Neil [11], 
Kakimoto et al.[3,4] and Mitchell [12]. The values of Davidson and O'Neil are 
on average 2.6 times larger than ours for nitrogen, but 0.88 of ours for air. 
The agreement of the present results with Kakimoto et al., whose experimental 
conditions are similar to the present ones, is rather good both in nitrogen and 
in air. The efficiencies of several bands calculated from the experiment of 
Mitchell using X-ray photons from 0.9 to 8.0 keV [12] are also listed. Though 
his excitation method and the target thickness are quite different from ours, 
it should be noted that the efficiencies at 800 hPa for air are in fairly good 
agreement with each other. 

The efficiency $° has been determined from the intercept at zero pressure of 
Eq.(2) in past measurements at low pressures [13] ^ 



*i(p) 



for air in ^ is plotted 



as a function of pressure in Fig. 6. The results of LS fitting to Eq.(2) are 
shown by solid lines. The determined $° from the intercept at zero pressure 
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Table 5 

Comparison of fluorescent efficiencies ^i{p) for N2 and air at p=800 hPa with 
those of D&OiDavidson and O'Neil [11], K:Kakimoto et al. cited from [3], [4] and 
M:Mitchell [12]. 



wave 

length 

nm 


Nitrogen 


Air 


Present 


D&O 


K 


M 


Present 


D&O 


K 


M 


xlO-4 


xlQ-^ 


337 


1.90 ±0.06 


5.20 


1.87 


2.45 


2.34 ±0.11 


2.10 


2.1 


3.38 


354 


0.093±0.010 


0.47 






0.284±0.047 


0.32 






358 


1.23 ±0.07 


3.70 


1.35 




1.73 ±0.17 


1.50 


2.2 




376 


0.182±0.012 


0.37 






0.489±0.074 


0.30 






381 


0.425±0.050 


1.40 




0.585 


0.58 ±0.10 


0.52 




0.79 


391 


0.086±0.010 


0.10 




0.114 


0.598±0.040 


0.70 


0.84 


0.64 


394 


0.037±0.015 








0.123±0.065 


0.05 






400 


0.079±0.007 


0.20 






0.249±0.037 


0.18 






406 


0.140±0.029 


0.43 




0.180 


0.226±0.017 


0.18 




0.227 


428 


0.028±0.006 


0.11 






0.220±0.039 


0.27 







and p' from the slope are in good agreement with the values listed in Table 
2. The dashed lines are from Mitchell [12]. The agreement is good for the 391 
and 406 nm bands, but not good at 337 nm and 381 nm. Open circles are 
from Davidson and O'Neil [11] and the asterisks are from Kakimoto et al. [4]. 
Their agreement with the present results for the 337, 354, 358, 381, 391 and 
428 nm bands are good, while there are slight differences for the 400 and 406 
nm bands. 

As shown in Fig. 6, the present results are in fairly good agreement with 
Davidson and O'Neil [11] and Mitchell [12]. The values oi p' and $° for the 
391 nm band from Mitchell are in good agreement with Hirsh et al. [13] mea- 
sured at low pressure. According to Mitchell, for the 391 nm band there is 
a significant deviation from Eq. (2) at pressures higher than 100 torr (133 
hPa). This deviation results from a three-body deactivation process and can 
be expressed as 



■^{1 + Kp + Dp^ 



:iii 



$(p) $0 

where $°=0.53%, K=1.08torr-i (0.810hPa-i) and D=4.4xl0-^torr-2 (2.48xlO-%Pa-2) 
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Fig. 6. The reciprocal of the efficiency ^^ in air is plotted as a function of pressure. 
The best-fit functions are shown by solid lines. The relations of Mitchell [12] are 
shown as dashed lines in some bands. Open circles are from Davidson and O'Neil 
[11]. Asterisks are from Kakimoto et al. [4] 

for the 391 nm band, which is shown by a dashed curve in the figure. There- 
fore we may conclude that our efficiencies at high pressures coincide with 
other data within experimental error, even if the present p' is different from 
that determined at low pressure by Hirsh et al. 

The relative intensity of each band in air is plotted in Fig. 7 with those from 
Bunner [2] and Ulrich et al.[14]. The incident energies of electrons are quite 
different (present:0.85MeV, Ulrich et al.:15keV Bunner:mainly 50keV), but 
the agreement may be within the experimental error of each experiment. 
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Fig. 7. Relative photon intensities in air at 1013 hPa in each band between SOOnni 
and 430 nm. The intensities are normahzed to 1.0 at 337 nm. 

We may, therefore, surely conclude that the present results can be applied to 
cosmic ray experiments, where photon yields at high pressures above 100 hPa 
are important. 



7 Remarks to the energy determination of UHECR 



7. 1 Effects to the UHECR experiments so far 



In order to estimate the effects of the present results to the energy determi- 
nation of the UHECR experiments so far made, let us compare the observed 
number of photons from the extensive air showers using the present photon 
yields with those used by the HiRes experiment [15]. The HiRes experiment 
currently uses a combination of results from Kakimoto et al. [4] and Bunner 
[2]. 

The conditions for the present calculation are similar to a previous one [16] 
and are summarized below. 

• Air shower simulation code used: CORSIKA 6.020 [17,18] with QGSJET 
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model [19]. 

■ Primary particle=proton 

■ Energy E=l0^^eV and lO^OeV 

■ Zenith angle ^ = and 60° (In the case of ^ = 60°, the shower axis is on 
the plane perpendicular to the line of sight.) 

■ In each combination of above parameters, 30 events are simulated to get 
the average longitudinal development curve. 

The number of emitted photons is calculated from the total energy deposit 
of shower particles, SdE'/dx, in each step for sampling of longitudinal devel- 
opment, using Eq. (8). The energy release provided by CORSIKA was used 
for this calculation, in which the contribution by particles below simulation 
energy threshold is carefully taken into account [20]. 
The observation height is m a.s.l. 

Emitted photons in each step width of depth are attenuated by Rayleigh 
scattering with the following transmission factor 



Tr = exp 



1 3^1 ^ ^2! /400[nm] 



Xn \ A 



:i2) 



where Xr=2974: g/cm^, and xi and X2 are the slant depths of two points. 
Photons are attenuated by Mie scattering with 



1m = exp 



Lm cos 9 






400|nm|\ ^^.^^ 



where scale height Hm = 1.2km and the horizontal attenuation length Lm = 
25km. hi and /12 are the heights of the emission point and the detection point 
of the light, respectively. 

• Total number of observed photons is calculated by adding photons from 
each step of depth. 

• US standard atmosphere 1976 [7] is used for the altitude dependences of 
density and temperature. 

• Wavelength dependences of the HiRes filter transmission and quantum ef- 
ficiency (QE) of the HiRes Photomultiplier tube (PMT) are taken into ac- 
count in some cases. 

In Fig. 8, the ratios Rn of the number of total observed photons with the HiRes 
photon yields {NniRes) to that with the present ones (Npresent) are plotted for 
various conditions as a function of distance (horizontal distance to the 2 km 
point a.s.l. along the shower track). NniRes is smaller than Npresent by —11% 
to —15% depending on the distance to the shower due to Rayleigh scattering 
(Fig. 8 (a)). In the case of the inclined shower, the change of R^ with distance 
is smaller than for the case of the vertical shower, and the difference is between 
— 11% and —13% (Fig.8(b)). Rn changes little if we include the effect of Mie 
scattering (with its weak wavelength dependence) (Fig.8(c)) and is almost 
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Fig. 8. Comparison of observed total photon number in various conditions. 
^N = jy""'^"" is plotted as a function of distance, (a) The effect of the primary 
energy for the showers with 9 = 0°. Circles with a solid line are for log£'(eV) = 20 
and triangles with a dashed line for log£'(eV) = 19. Only Rayleigh scattering is 
taken into account, (b) The effect of zenith angle for log£^ = 20. Only Rayleigh 
scattering is taken into account. A solid line is for 0° and a dashed line for 60°. (c) 
The effect of Mie scattering for log E = 20 and ^ = 0° . Circles are for the case that 
only Rayleigh scattering is taken into account, while triangles are for the case that 
Mie scattering is also considered, (d) The case of HiRes filter transmission and PMT 
QE for the showers with log-E = 20,^ = 0°. Both Rayleigh and Mie scatterings are 
taken into account. Circles are the same with those in (c). 

independent of the primary energy (Fig. 8 (a)). However, if the HiRes filter 
transmission and the QE of the PMT are taken into account, R^ becomes 
closer to unity, —4% to —7%, depending on the distance (Fig. 8(d)). This is 
because the transmission coefficient of the broad band filter used in HiRes 
drops below 320nm and above 380nm where the difference of photon yields 
between the present experiment and those assumed by HiRes are relatively 
large. 

We need more a detailed analysis based on the experimental conditions to 
infer the individual energy estimation of cosmic rays, taking into account 
the wavelength dependent items (photon yields, scattering effects, detector 
performance, etc.). 
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1.2 Application of the present photon yields 



In UHECR fluorescence experiments, the primary energy is estimated from 
tlie calorimetric energy, Ecai, witli a correction for missing energy, E^iss, that 
carried by neutrinos and muons, and that lost due to nuclear excitation [21]. 
Ecai may be determined from the experiment from the path-length integral 
multiplied by the mean ionization loss rate, a, over the entire shower as 



E,ai = a J N,h{X)dX , (14) 



where Nch{X) is the number of charged particles in the shower as a function 
of atmospheric depth X in g/cm^. Song et al. [21] showed that this technique 
provides a good estimate of the primary energy of cosmic rays with q;=2.19 
MeV/(g/cm^) by using the CORSIKA air shower simulation program. 

Their arguments would be accepted as far as Nch{X) can be determined un- 
ambiguously. In practice, however, it is not straightforward to convert the 
observed number of photons for each angular bin in the camera of the detec- 
tor to Nch{X). This is because the photon yield depends on the temperature 
and density of air along the trajectory of the electrons, and most particles of 
low energies are not traveling parallel to the shower axis. Alvarez-Muhiz et 
al. [22] studied the ratio of the average track length traveled by the shower 
particles in a some depth interval and that projected onto the shower axis. 
They showed that the ratio depends on the shower age and is 1.18 at shower 
maximum. That is, E^ai is possibly overestimated, if Eq. (14) is used without 
path length correction. 

Instead of estimating Nch{X), Dawson [15] has proposed to use the energy 
deposited in the atmosphere by the shower per g/cm^ of depth. This is de- 
termined from the number of photons, L[X), after correcting the attenuation 
of photons due to Rayleigh and Mie scattering and subtracting Cherenkov 
contamination. That is, the energy deposited by the shower in the grammage 
interval AX, AEdep{X), is expressed by 



AE,,,iX) _ 1 ( L,iX)ihu), \ 

AX AX^[ MX) J' ^ ' 

where Li{X) is the fraction of the flux L{X) in different wavelength bin i and 
(/iz/)j is the photon energy of bin i. 

Then the calorimetric energy is estimated from 
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E..= l^^^dX, (16) 



It would be much more realistic to use Eqs. (15) and (16) to estimate the 
primary energy. $i(X) is a function of pressure {=pRAirT) at depth X and is 
expressed by Eq. (2). Recommended values of $° and p' are listed in Table 2 
for each wavelength. 



8 Conclusion 



Photon yields have been measured in fifteen wave bands as a function of 
pressure, for nitrogen and dry air excited by electrons of an average energy 
of 0.85 MeV. The pressure dependencies of 15 wave bands between 300 nm 
and 430 nm have been determined with our own measurements. The total 
photon yields between 300 nm and 430 nm are 22.20±0.56 and 4.05±0.14 per 
meter per electron at 1013 hPa and 20°C for nitrogen and air, respectively. 
If we restrict the wave bands up to 406 nm, the corresponding values are 
21.69±0.55 and 3.81±0.13. The systematic error in the measurement is 13 %. 

From the pressure dependence of photon yields, their temperature and den- 
sity dependencies of each band are determined for application to the energy 
estimation of UHECRs by the fluorescence method. It would be much more 
realistic to use energy deposition rather than the number of charged particles 
for each angular bin in the camera of the detector to estimate the primary 
energy. <I>° and p' in each wavelength are given in Table 2 for application to 
estimate the energy deposition. 
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